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A growing number of shock compression experiments, especially those involving laser compression, 
are taking advantage of in situ x-ray diffraction as a tool to interrogate structure and microstructure 
evolution. Although these experiments are becoming increasingly sophisticated, there has been little 
work on exploiting the textured nature of polycrystalline targets to gain information on sample 
response. Here, we describe how to generate simulated x-ray diffraction patterns from materials 
with an arbitrary texture function subject to a general deformation gradient. We will present 
simulations of Debye-Scherrer x-ray diffraction from highly textured polycrystalline targets that 
have been subjected to uniaxial compression, as may occur under planar shock conditions. In 
particular, we study samples with a fibre texture, and find that the azimuthal dependence of the 
diffraction patterns contains information that, in principle, affords discrimination between a number 
of similar shock-deformation mechanisms. For certain cases we compare our method with results 
obtained by taking the Fourier Transform of the atomic positions calculated by classical molecular 
dynamics simulations. Illustrative results are presented for the shock-induced a—e phase transition 
in iron, the a-ui transition in titanium and deformation due to twinning in tantalum that is initially 
preferentially textured along [001] and [Oil]. The simulations are relevant to experiments that can 
now be performed using 4th generation light sources, where single-shot x-ray diffraction patterns 
from crystals compressed via laser-ablation can be obtained on timescales shorter than a phonon 
period. 


I. INTRODUCTION 

The response of matter to shock compression has been a subject of study for more than a half a century 1-0 , 
and early in this field of research it was recognised that the high transient stresses to which materials could be 
subject caused them to yield and deform plastically, and, in certain circumstances, to undergo rapid polymorphic 
phase transitions'. An understanding of the pertinent physics operating at the lattice level that underpins the 
material response has long been sought, and has been a strong motivator in the development of lattice-level and 
structural measurement techniques with sufficient temporal resolution to interrogate detailed material response during 
the deformation process itself. Of particular relevance to the work presented here are recent advances in ultra-fast 
x-ray diffraction. The first diffraction patterns of crystals subjected to shock compression taken during the passage 
of the shock itself had exposure times of ten of nanoseconds 11 Ml . Since this pioneering work significant progress has 
been made by the use of x-ray sources based on diode-technology 1 L synchrotron emission 1 and the use of x-rays 
emitted from plasmas created by irradiation with high-power lasers 1 ' 1 ’ 1 . More recently however, the development 
of so-called 4th generation light sources, such as the Linac Coherent Light Source (LCLS), has allowed high-quality 
single-shot diffraction patterns to be obtained in just a few tens of femtoseconds - freezing motion on a timescale 
faster than the shortest phonon period in the system. 

The transient x-ray sources mentioned above have been used to further our understanding of the response of 
both single-crystals and polycrystalline matter to both shock and quasi-isentropic compression 1 ’ 1 ’ ■ “ . A num¬ 
ber of x-ray diffraction techniques have been developed to monitor material response, including divergent beam 
geometry 50 , white-light Laue diffraction l , Debye-Scherrer diffraction “ 1 and the use of energy-dispersive single¬ 
photon counting 32 . While studies of both single-crystals and polycrystalline matter have been undertaken, in the case 
of polycrystalline samples very little attention has yet been paid to the role of texture in these uniaxial compression 
experiments - that is to say the distribution function of grain orientations within a particular polycrystalline specimen. 
A number of manufacturing methods, such as rolling and epitaxial growth, result in characteristic textures due to the 
way in which the materials have been processed. This texture has an influence on a range of physical properties such as 
strength, electrical conductivity and wave propagation . Therefore texture plays an important role in understanding 
material response. 
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Furthermore, whilst the texture of a sample will influence its response to rapid uniaxial compression, it will also 
have a profound influence on the way in which the sample diffracts. As certain grain orientations are more likely to 
occur than others, the intensity of a particular Debye-Scherrer ring (corresponding to a certain set of Miller indices) 
will have a strong azimuthal dependence, and this dependence can in turn be used to extract texture information. 
Indeed, significant static studies of the texture of polycrystalline samples have been undertaken with synchrotron 
sources for many years. 3,1 ' . Wenk and co-workers provide an overview of the use of synchrotrons in such texture 
analysis ’ ’. While the texture of a material is often represented by a set of pole figures which can be measured 
directly via x-ray diffraction, prediction of anisotropic material properties requires knowledge of the full orientational 
distribution function (ODF). The ODF gives the probability of a crystallite having an given orientation, therefore 
provides a complete description about the texture of the sample. Pole figures, being a 2D projections of the ODF, 
result in some texture information being lost, although methods have been developed to obtain an approximate ODF 
from pole figures 3 ’. 

Given that the preferred orientation defined by texture links both the diffraction patterns observed, and the sample 
response, it is logical to question whether specific information can be gleaned via in situ diffraction studies of samples 
with known, well-defined texture. For example, bulk rotations of the crystal lattice, or changes in the crystal structure, 
such as Martensitic phase changes, will result in an altered texture that could be used to distinguish between different 
mechanisms of atomic rearrangement. This reorientation has been observed in previous work using both neutron 
sources 111 and synchrotrons 3 K ', although only at relatively modest pressures compared with those we are interested 
in here. 

Ideally one would wish to understand the detailed response to shock compression of samples as a function of their 
ODF, and to then predict the resultant diffraction patterns. We outline in section II the method by which this could 
in principle be done. However, in terms of the results of particular simulations, our present goal is more modest: we 
restrict ourselves to crystals that are highly fibre-textured - that is to say all of the grains within the sample initially 
have very similar orientations with respect to a single axis (the fibre axis), only deviating slightly in angle from that of 
a high-symmetry direction. As all of the grains have similar orientations along a given axis, which we also take to be 
the axis of compression, we might expect that we can, to a reasonable approximation, model certain aspects of their 
response to shock or quasi-isentropic compression using single-crystal parameters. In particular, when determining 
how the material will respond to compression and shear, we assume that the elastic stresses, supported by elastic 
strains, can be calculated by using elastic constants appropriately chosen to mimic single crystal response along the 
pertinent directions. In terms of x-ray diffraction, however, the finite range of orientations determined by the ODF 
is such that a monochromated, non-divergent incident x-ray beam can diffract from a reasonably large subset of the 
grains, both in the shocked and un-shocked case. A similar approach to that outlined above has recently been used 
to observe, via femtosecond diffraction, the ultimate compressive strength of copper subjected to shock compression 
at strain-rates of order 10 9 s -113 . We shall show that breaking the symmetry of the system, such that the direction 
of the incident x-rays and the target normal (parallel to the compression direction) are no longer antiparallel, allows 
us to determine specific information about the system under study. 

The layout of the paper is as follows. In section II we describe the method by which we determine the analytic 
diffraction pattern as a function of compression, and the angle that the sample normal makes to the x-ray beam. 
The background to the molecular dynamics simulations are outlined in section III before, in section IV providing the 
results for simulations for the a-e phase transition in iron, the a-co transition in titanium, and deformation due to 
twinning in tantalum. We then conclude with a discussion of the results, and suggestions for experiments and further 
work. 


II. METHOD 

We have previously shown how to calculate the position of Debye-Scherrer rings from a polycrystal strained by an 
arbitrary deformation gradient in the Voigt limit 11 , and for the sake of completeness we restate briefly the results here. 
We envisage a planar sample orientated such that its surface normal n lies along the z' direction. As mentioned above, 
we allow symmetry to be broken, such that the direction of incidence of x-rays for wavevector k 0 can be non-parallel 
to n, and along the direction z. The rotation matrix R transforms from the unprimed (x-ray) to primed (target) 
co-ordinate system. Debye-Scherrer diffraction occurs when the Laue condition is met: G = R 7 G' = k — ko, where 
k is the wave-vector of the diffracted x-ray. Deformation of the crystal in real space is represented by the arbitrary 
deformation gradient F applied in the target (primed) co-ordinate system, which corresponds to the analogue in 
reciprocal space of (F 7 ) _1 . With these relations we see that the reciprocal lattice vector for an unstrained crystal in 
x-ray co-ordinates, Go, is transformed under deformation to a new reciprocal lattice vector G, where 

G 0 = R t F t RG = aG, 


(1) 
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where a = R 2 F 2 R. By noting that if a plane is to diffract once the sample has been strained, it must meet the Laue 
condition, and that for the unstrained case, |G 0 | = it is possible to solve equation 1 to yield the Bragg angle as a 
function of azimuthal angle around the z-axis, <j>, for an arbitrary deformation gradient - that is to say we know the 
direction of the scattered wavevector, k, and hence can determine via simple ray-tracing where the diffraction will 
impinge on a distant detector. In Higginbotham et al. n , diffraction patterns for numerous deformation gradients and 
target geometries (with respect to the x-ray beam) are shown. However, in these cases the sample was assumed to be 
isotropic in texture, and thus satisfying equation 1 was deemed a sufficient condition for diffraction to occur. 

However, a non-isotropic ODF places further constraints on the possibility of diffraction for a given Bragg angle 
and azimuthal position around the Debye-Scherrer ring, as the ODF provides the measure of the probability of finding 
a crystal with a given G 0 existing in the original unstrained sample. The route forward for simulating diffraction 
for a crystal with known original ODF under a known arbitrary deformation gradient is now clear: we use the ODF 
to determine the probability of a given crystallite having the appropriate Go existing within the sample, we then 
determine G from equation 1 (and hence k from the Laue condition), and use ray-tracing to propagate the diffracted 
beam to the detector assuming that the intensity is proportional to the probability of finding the original Go, as 
determined by the ODF, and taking into account factors such as multiplicity. 

Although the above approach is general, within the rest of this paper we restrict our study to the particular case 
of a simple fibre texture where the crystallites have nearly identical orientation in the axial direction, but close to 
random radial orientation. Our motivations for this are due to the fact that this allows us to treat the mechanical 
response of the polycrystal to be well-approximated by that of a single crystal with orientation aligned with the fibre 
axis and that the technique of fibre diffraction under ambient conditions is well established 2 ’ . Furthermore, recent 
experiments using femtosecond x-rays created with 4 th generation light sources to diffract from uniaxially compressed 
samples have employed targets with this type of texture, 18 and the thin films that have hitherto been used in these 
experiments often grow with such preferential orientation. 

Fibre textured samples have a greatly simplified ODF. For the case of perfect fibre texture, each crystallite has 
a single crystallographic direction, the fibre orientation, vj., associated with the reciprocal lattice vector (hi,ki,li), 
which for all crystallites are aligned parallel to the sample normal, n. Each grain is then deemed to possess a random 
orientation when rotated about the axis, vj. If we consider a particular plane within a crystallite, with miller indices 
(h 2 , ^ 2 , h) (which are the same set of miller indices associated with Go) to which the reciprocal lattice vector is V 2 , 
then the value of vy • v 2 will be a constant. However, when z and z' are non-parallel, vy • G 0 varies as a function 
of azimuthal angle around z. Thus, for perfect fibre texture the Debye-Scherrer pattern is not a ring pattern, but an 
array of points defined by simultaneously satisfying equation 1, as well as the condition vy • Go = Vi • v 2 . Note that 
this condition does not take into account multiplicity, and therefore one must consider this condition for each member 
of v 2 in the {/i 2 , ^ 2 ,^ 2 } family, since they do not necessarily result in the same values of Vi • v 2 . 

However, real fibre-texture samples contain crystallites that are not perfectly aligned axially, and the volume fraction 
of crystallites with reciprocal lattice vectors Vi, P(vi), will be a rapidly decreasing function of vy ■ v n , where now 
v n =< vi >. We further assume that for a given volume fraction with a certain vj, the directions v 2 of the normals 
to the planes with miller indices [^. 2 ,^ 2 , ^ 2 ] within the crystallites will be arranged randomly azimuthally around the 
axis Vi, such that the possible orientations of v 2 are simply constrained by v 2 ■ vi being equal to the value it would 
take for a single crystal: that is to say P(v 2 ) oc P(vi) subject to the v 2 • Vi constraint. The diffraction condition is 
once more defined by simultaneously satisfying equation 1, as well as the condition vy • Go = Vi • v 2 ., but now the 
intensity of the diffraction is proportional to P(vi). As the possible directions of vi can now vary over a constrained 
range, our pattern is a series of arcs, rather than points. For the sake of simplicity in our simulations here we assume 
that the volume fractions of crystallite orientiations are uniform over a small range of angles, i.e. P(v 1 ) = C, a 
constant, for | arccos(vi • Go) — arccos(vi ■ v 2 )| < S , and P(vi) = 0 for | arccos(vi ■ Go) — arccos(vi • v 2 )| > S. 
This method effectively finds the intersection between the Ewald and Polanyi spheres 42-44 , however for the case of 
anisotropic strains, the Polanyi sphere is deformed into an ellipsoid. Note that by changing the x-ray energy or the 
sample orientation, which varies the size of Ewald sphere or rotates the Polanyi sphere respectively, the position of 
the arcs on the Debye-Scherrer ring also change, and that by varying these parameters, different parts of reciprocal 
space can be interrogated. 


III. MOLECULAR DYNAMICS SIMULATIONS 

One advantage of the fibre texture discussed above is that one knows (to within the texture width) the crystallo¬ 
graphic orientation of grains with respect to an applied planar compression front. This is particularly important as 
crystallite orientation can drastically alter the material response under uniaxial loading 2 ' ’ - .In the case that grain 
size is comparable to sample thickness , one can approximate the response of the sample as being close to that of a 
suitably oriented single crystal. This is particularly pertinent if one wishes to compare results with those of molecular 
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FIG. 1: (Color online) A simulated ray trace of diffraction from a [001] fibre textured polycrystalline HCP Fe formed under 
shock compression. The blue overlaid lines show the positions of arcs from uniaxially compressed BCC, while the red lines 
show the positions of arcs from HCP with OR [001] bcc || [2ll0] hcp , [110] bcc || [0002] hcp . For clarity the x-ray energies used to 
trace the bcc and hep overlays were offset by 1%. The corners of the detector are located at a scattering angle 26 = 74.2°. 


dynamics, where state of the art polycrystalline simulations are still generally limited to grain sizes of 5-100nm 
far below the grain sizes of typical experimental samples. 

In order to relax the requirements on computational power, we present a method of simulating the response of a fibre 
textured target by manipulation of a single crystal simulation. We do this by first performing a 3D Fourier Transform 
(FT) of the computed electron density of the single crystal . This provides us with a momentum space representation 
of the lattice which describes the allowed scattering vectors for diffraction. 

Working with the intensity of this FT, we first note that any polar dependence around the compression axis can 
be neglected due the random rotational distribution of grains in a fibre textured sample. Considering a cylindrical 
geometry, we therefore produce a 2D representation of the FT, which flattens the data into its (/?, z) components, 
effectively integrating around <f>. In the case of a perfect (zero texture width) fibre texture, this 2D representation 
correctly describes all scattering. 

For the case of finite texture width, one can imagine that the misorientations of the grains are simply related to 
a rotation about the origin of the 2D representation, and so to mimic the width we sum rotated representations 
for angles between ±5. This new representation necessarily still retains the cylindrical symmetry required for fibre 
texture, but correctly accounts for the distribution of grain alignments. Of course, as S becomes larger than a few 
degrees, the underlying assumption that all grains react in a similar manner to a well aligned single crystal will break 
down. However, for this paper we will assume suitably small texture widths of ss 5° where this approach works well. 
One can now raytrace simulated diffraction patterns directly from this 2D representation by only considering the 
( p , z ) component of the scattering vectors expressed in this cylindrical target geometry. 


IV. RESULTS 

We now present results for three different fibre-textured polycrystals subject to uniaxial compression: iron with 
a [001] fibre texture, titanium with [0001] fibre-texture, and tantalum with [001] and [011] fibre-textures. For all 
three crystals we make assumptions about the deformation mechanisms, and use the approach of section II to predict 
diffraction patterns. For two of the cases - iron and tantalum - we also compare our calculations with diffraction 
patterns simulated using the method of molecular dynamics, as outlined in section III. 
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A. a—e Phase Transition in [001] Iron 

The a-e phase transition in iron is an example Martensitic transformation, characterised by a collective movement 
of atoms across distances that are typically smaller than a nearest-neighbour spacing. These type of transitions 
are well suited to laser compression studies, since the timescales on which they occur are comparable to the short 
pulses that can be attained in laser experiments. Importantly, for these non-diffusional transitions, an orientational 
relationship (OR) exists between the two phases. While the OR does not uniquely determine the mechanism by which 
the phase transition occurs, it can significantly reduce the number of candidate mechanisms. The ability to measure 
the OR in situ is therefore highly desirable. 

As an example of determination of an OR, we take the case of [001] oriented iron, where the phase transition 
mechanism is well understood. Molecular dynamic simulations undertaken by Kadau et al. aimed to understand 
iron’s bcc-hcp shock induced phase transition 1 . The results of these investigations were later reproduced with 
remarkable fidelity in experimental x-ray diffraction studies ' . In both MD and experiment, the OR was described 
as [001] bcc || [2110] h , [110] bcc || [0002] hcp " 1 . One can therefore consider this OR as reorienting a fibre textured sample 
from (002) bcc to (2ll0) hcp . 

Following Kadau, we simulate a 100x100x800 cell (288x288x230lA) iron single crystal shocked along the [001] 
direction by 0.7km s -1 piston using the same Voter-Chen potential used in Kadau’s work . For this piston velocity, 
the material does not reach the 18.4% uniaxial compression needed to create ideal HCP, instead reaching only 
13.8% 25 , resulting in an anisotropically strained HCP structure. A 3D FT was performed on a section of the material 
behind the shock front. The FT was modified to mimic that of a fibre textured polycrystal, using the method 
described in section III. Figure 1 shows the resultant ray trace 1 for a detector in transmission geometry, using 
a 12 keV x-ray source and with the sample normal rotated at an angle 30° to the incoming x-rays. The overlaid 
red lines show the predicted diffraction pattern (using the methods of section II) of strained HCP iron described 
above, with a c/a ratio of 1.73, a texture direction along [1120] h and a textured width of 5°. The blue lines 
show the pattern from 13.8% uniaxially compressed BCC. As expected, there is agreement between the raytrace 
and predictions from the molecular dynamics simulations, supporting the validity of the approach outlined in section II. 



FIG. 2: (Color online) The predicted diffraction pattern using 7.5 keV x-rays from the uj phase of a shocked [0001] fibre textured 
Ti sample, with the sample normal rotated at an angle 25° to the incoming x-rays. The green arcs labelled A and B correspond 
to the {1011} and {1120} sets of planes respectively. The corners of the diagram are located at a scattering angle 29 = 70.5°. 
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TABLE I: The three ORs that correspond to various proposed mechanisms of the a — cu phase transition in [0001] Titanium 

Variant Orientational Relationship Ref 

i (oooi)ji(ioiiL , [ioio]J|[ioiil 55 “ 57 

II (0001)„ II (1210)„ , [1210] q II [0001]„ 55 ’ 58 

Zong (0001) a || (10I0)„ , [10I0] Q || [1123]„ 47 


B. a—LO Phase Transition in [0001] Titanium 

The group IV transition metals titanium (Ti), zirconium (Zr) and hafnium (Hf) have a hexagonal close packed 
structure (a) under ambient conditions, but exhibit a Martensitic phase transition to another hexagonal structure 
(ui) under high pressure. Although the a-u transition is well established, the mechanism by which it occurs, and 
therefore the OR between the phases, is still not fully understood, especially under different loading conditions. A 
summary of the ORs for various proposed mechanisms is given in table I. The first two ORs were proposed by Usikov 
and Zilbershtein in a TEM study of statically compressed Zr and Ti, by arguing that the transition occurs via an 
intermediate /3 phase, and are usually referred to as Variant I and II. Earlier work by Silcock ' on the oj phase in Ti 
alloys proposed a different direct mechanism which corresponds to the Variant II OR. Later computational work by 
Trinkle demonstrated that a new mechanism, known as TAO-1, had a lower energy barrier than that proposed by 
Silcock. This new mechanism produced the Variant I OR. Experimental work by Song and Gray 5 ' observed the OR 
(0001) a || (1011)^, [1010] a || [1123]^, although independent analysis re-examining this data led to the conclusions that 
is actually a subset of Variant I . We have therefore associated Song and Gray’s work with Variant I in table I. 

Molecular dynamics simulations performed by Zong et al. found that Ti shocked along the [0001] direction resulted 
in a mechanism that gave an OR differing from both Variant I and II, as listed in the table. Note that each OR results 
in a different fibre texture direction for the w-phase for uniaxial compression of an initially fibre-textured a-phase 
target, and thus a target with the normal tilted with respect to the incident x-ray beam is well suited to provide some 
discrimination between variants, therefore providing some insight into possible transformation mechanisms. 

The work by Zong 1, found that Ti shocked along the [0001] direction by a piston with a velocity of 0.9 km s _1 
resulted in transformation to the ui phase, with lattice parameters = 4.61 A and = 2.82 A. Using these values, 
we calculate the diffraction pattern from the ui phase of a shock-compressed [0001] sample of Ti with an angular 
texture width of ±5°. The sample normal is set at angle 25° to the incoming x-rays, which have an energy of at 
7.5 keV. The predicted diffraction patterns for each of the possible variants are shown in Figure 2. The blue, red and 
green lines correspond to the diffraction from u material of Variant I, Variant II and the Zong OR respectively, while 
the dotted black line corresponds to the diffraction from an untextured polycrystalline sample. For clarity, the blue 
lines have been slightly offset outside the true Debye-Scherrer rings, while the red lines have been slightly offset inside. 
A clear difference can be seen in the diffraction patterns for the different variants, thus allowing the OR and hence a 
subset of mechanisms to be determined by the azimuthal position of the diffracting arcs around the Debye-Scherrer 
ring. It is important to note that the ability to discriminate clearly between all three ORs is not guaranteed, and 
relies on a judicious choice of both x-ray energy and tilt angle. 

Figure 2 also demonstrates that lines with similar d-spacings do not necessarily appear at similar azimuthal positions. 
For example in an ideal cj crystal, the {1011} and {1120} planes have the same spacing and are therefore completely 
unresolvable by powder diffraction from an untextured sample. However, since these planes form different angles to 
the sample normal, within a textured sample their corresponding arcs appear at different azimuthal angles around 
the Debye-Scherrer ring, allowing them to be resolved. This is shown in Figure 2, where the green arcs labelled A 
correspond to the {1011} set of planes, while the green arc labelled B corresponds to the {1120} set of planes. By 
resolving these two arcs, we are able to gain information that cannot be obtained via powder diffraction from an 
untextured sample; in this case on any small departure from the ideal c/a ratio. We note that it is only possible to 
resolve lines with similar d-spacings if the angle between G-vectors of each of the planes and the sample normal is 
significantly different. 


C. Twinning in [001] and [Oil] Tantalum 

Tantalum provides an appealing case to study, owing to its multitude of competing plasticity mechanisms (a 
combination of dislocation and deformation twinning mediated responses). This is of particular interest, as Debye- 
Scherrer diffraction in polycrystalline samples with completely random texture cannot distinguish between slip and 
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FIG. 3: (Color online) Simulated diffraction from a [001] fibre textured polycrystalline Ta under shock compression. Overlaid 
are the predicted diffraction patterns for: (a), the untwinned case, with elastic strains of tt = —0.029, ez = —0.136 shown in 
green and (b) the twinned case, with elastic strains of et = —0.052, ez = —0.109 , shown in red. The corners of the detector are 
located at a scattering angle 26 = 66.3°. 


twinning. Most experimental work into twinning under uniaxial dynamic compression has been performed using 
either explosive lenses or gas guns' . However, to date, time resolved laser diffraction experiments have failed to 
yield any evidence of twinning in situ , although residual twinning has been observed in laser driven shock recovery 
experiments' 11,65 . 

As with phase transitions, the reorientation of the lattice caused by twinning will result in a different crystallographic 
direction being oriented along the fibre direction. This may, in turn, lead to a signature in the diffraction pattern. A 
similar technique using neutron diffraction has been used to observe twinning in magnesium, which occurs at much 
lower pressures than in tantalum . 

Molecular dynamics simulations performed by Higginbotham et al. have predicted a significant twinning fraction 
in [001] tantalum under shock compression . In that work, the sample was found to be almost completely twinned 
when compressed by a piston with a velocity of 0.9 km s — , corresponding to a uniaxial compression in the elastic 
wave by 18%. The authors noted that after an initial uniaxial compression of 18.4%, twinning could be achieved by 
shuffling alternating (112) planes in the (111) direction. They therefore proposed this to be the mechanism by which 
the twinning occurred, with the material reaching its final state via elastic relaxation, although they caution that, 
given the relatively simple nature of the potential used, they do not claim to exactly model what will occur in practice 
in shocked Ta. However, the observed shuffling provides a possible mechanism for how twinning of bcc materials may 
occur under shock compression. 

We repeated the simulations in this work, using a 100x100x800 cell (330x330x2644A) Ta single crystal, modelled 
using the EFS potential", and shocked along the [001] direction by a piston travelling at 0.9km s -1 . The per atom 
structure factor (PASF)" was used to distinguish between twinned and untwinned material in the plastically deformed 
material behind the shock front. A 3D FT was performed on a stable region within the plastically deformed material 
behind the shock front, as well as the twinned region separately, and both were modified in the way described in section 
III. The diffraction pattern was simulated assuming a 12keV x-ray source, with the angle between the incoming x- 
rays and the sample normal being 25°. The results are shown in Figure 3. For the case of slip, while small rotations 
of the crystal lattice have been observed '", no large reorientation is expected, and thus the fibre orientation will 
remain close to [001]. However, by comparing the observed and predicted positions of the diffraction arcs for a [001] 
textured sample including the strains described above (shown in Figure 3a), it is clear that there has been significant 
reorientation of crystallites within the sample, indicative of twinning. 

In Figure 3b we plot the predicted diffraction pattern using the methods of section II assuming that the sample 









FIG. 4: (Color online) Simulated diffraction from a [Oil] fibre textured polycrystalline Ta under shock compression. Overlaid 
are the predicted diffraction patterns for: (a), the untwinned case, with elastic strains of tt = — 0.045, e; = —0.075, shown in 
green, (b) the twinned case, with the [922] direction close to the compression axis, shown in red, and (c) both the untwinned 
and twinned case together. The corners of the detector are located at a scattering angle 26 = 51.3°. 


has undergone twinning, and is subject to the longitudinal and transverse strains noted above. In order to find the 
diffraction pattern resulting from the shuffling mechanism above, one must consider how this affects the lattice vectors 
of the crystal. The sample is first compressed uniaxially along the z axis by 18.4%. The shuffling then has the effect 
of causing the crystal to be reflected in the (112) twinning plane. Note that as crystal is uniaxially compressed, this 
causes the [001] direction to be reflected to the [111] direction of the compressed crystal, rather than the [221] direction 
expected under hydrostatic conditions (see Figure 6 of Higginbotham et al. ’). Lastly, the twinned material returns 
towards the hydrostat, by relaxing along the longitudinal direction and compression along the transverse directions 
in the new rotated coordinates. The final longitudinal and transverse strains were measured to be -0.109 and -0.052 
respectively. However, in the textured sample, the effect of a finite texture width must also be considered. In this case, 
the initial compression occurs at slight angle to the texture direction, which results in a slightly different orientation 
of the twinning plane. The crystal is then reflected in this altered twinning plane, and relaxes as before. To create 
the predicted diffraction pattern, the lattice vectors of many crystallite orientations within the desired texture width 
were deformed by the method given above. These were then used to calculate deformed reciprocal lattice vectors, and 
thus the resultant diffraction pattern. 

Excellent agreement can be seen between the analytic solution, and the MD simulation, demonstrating that twinning 
has occurred, although slight differences can be observed which are due to the small angle assumptions used in section 
III to simulate the 3D FT of a fibre textured target. Additionally, there are some very weak arcs in the data 
corresponding to plastically compressed material (Figure 3a). The ratio of intensities of lines from twin and slip 
deformed material is indicative of the twin fraction. 

While most theoretical work on twinning in Ta has concentrated on the [001] direction, recent MD studies by Ravelo 
et al. have suggested shocking along the [011] direction may be more favourable for deformation twinning, due to 
the lower observed shear stress threshold for twin nucleation . This agrees with gas gun recovery work on Ta single 
crystals, which found a significantly higher twin volume fraction in shocks along [Oil], compared to the [001] and 
[111] directions' . Furthermore, as the [011] direction is the preferred direction for epitaxial growth of fibre textured 
thin films, this direction is particularly well suited for this technique. The [Oil] orientation exhibits two types of 
{112} (111) twin systems, which result in different fibre orientations. Under hydrostatic conditions, the first type 
causes no change in the fibre texture, while the second causes a reorientation to the [411] direction. However, only 
the second type has non-zero Schmid factors, and it is therefore assumed that only twins of this type occur. It follows 
that diffraction arcs corresponding to a [001] fibre orientation are from untwinned material, while arcs corresponding 
to fibre orientations close to the [411] type directions are from the two twinned variants. 

We repeated the simulations in this work, using a 100x100x800 cell (330x330x2644A) Ta single crystal, single 
crystal, modelled using Ravelo’s Tal EAM potential" 9 , and shocked along the [Oil] direction by a piston travelling at 
0.62 km s -1 . Again the PASF was used to distinguish between twinned and untwinned material, and a 3D FT was 
performed on each region separately. In the twinned material, the [922] direction of the compressed crystal is close 
to the compression axis, which is consistent with twinning after an initial uniaxial compression, similar to the [001] 
case. The diffraction pattern was simulated assuming a 10 keV x-ray source, with the angle between the incoming 
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x-rays and the sample normal being 45°. The results are shown in Figure 4. Since the twinning mechanism in [Oil] 
Ta is not well understood, the predicted pattern was produced for the structure found measured with the FT. Again, 
good agreement is seen between the analytic solution and the MD simulation, although in this case, there are strong 
arcs corresponding to both twinned and untwinned material, suggesting a significant amount of both are present in 
the sample. 


V. DISCUSSION 

The examples we have given above demonstrate that x-ray diffraction from uniaxially compressed fibre-textured 
targets can in principle yield information on deformation mechanisms, be they due to phase transformations or 
twinning. The breaking of the symmetry of the problem, by tilting the normal of the polycrystalline target with 
respect to the incident x-ray beam allows the encoding of such information in the azimuthal distribution of intensity 
in the Debye-Scherrer rings. We envisage that the methods that we have outlined here will aid in the design of 
experiments that have as their goal the elucidation of such mechanisms. It is worth considering, however, that the 
choice of initial fibre direction is important in determining what structural information can be extracted. In particular 
it should be noted that for [001] fibre oriented Fe and Ta samples, these orientation do not have the lowest surface 
energy, and thus are not the typical orientations in which thin polycrystalline foils of these materials grow. It may 
thus be that some effort is required to fabricate suitable samples. This is not an issue for the case of [0001] Ti or 
[011] Ta, which are usually grown with these textures. Beyond the four demonstration cases given above, it is clear 
that further work could concentrate on a variety of different samples and deformation mechanisms. In addition, we 
believe that the technique may have other advantages for the study of samples subject to shock or quasi-isentropic 
compression. Owing to the high strain rates present in such experiments, high dislocation densities or small grain 
sizes under phase transformation may ensue” ’ , resulting in broad diffraction peaks that are hard to resolve simply in 
terms of scattering angle, and thus would not necessarily be easily amenable to study by techniques such as Rietveld 
refinement. However, tilting of a target and separation of diffraction peaks azimuthally offers a possible route to 
finding structural solutions under the extreme pressures that can be obtained via laser-ablation. We believe that 
the technique we have outlined could have application to laser based shock and compression experiments that make 
use of emergent 4 th generation light sources, which offer incredibly bright, narrow bandwidth x-ray sources, with 
unprecedented temporal resolution. 
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